N 7 3- 12 8^1 . 

Reports of the Department of Geodetic Science 


Report No. 185 


ON THE GEOMETRIC ANALYSIS 
AND ADJUSTMENT 

OF OPTICAL SATELLITE OBSERVATIONS 


by 

Emmanuel Tsimis 


case mum 


Prepared for 

National Aeronautics and Space Administration 
Washington, D.C. 


Pi i ® ; 

V W'4; -- '■ •‘v ’ 1 

' 1. . 'v 'S'it 'Kl -I : ' ‘ ‘ ' 




•IV 


"* V. ■ ■ ■' 


,;.W, 


v>-v 


Contract No. NGR 36-008-093 
OSURF Project No. 2514 


■ f . 

. S’ ■ ' 



The Ohio State University 
Research Foundation 
Columbus, Ohio 43212 


August, 1972 


Reports of the Department of Geodetic Science 
Report No. 185 


ON THE GEOMETRIC ANALYSIS AND ADJUSTMENT 

OF OPTICAL SATELLITE OBSERVATIONS 

by 

Emmanuel Tsimis 


Prepared for 

National Aeronautics and Space Administration 
Washington, D. C. 

Contract No. NGR 36-008-093 
OSURF Project No. 2514 


The Ohio State University 
Research Foundation 
Columbus, Ohio 43212 

August, 1972 



PREFACE 

This project is under the supervision of Ivan I. Mueller, Professor 
of the Department of Geodetic Science at The Ohio State University, and 
it is under the technical direction of Benjamin Milwitzky, Project Manager, 
Geodetic Satellite Program, NASA Headquarters, Washington, D. C. The 
contract is administered by the Office of University Affairs, NASA, 
Washington, D. C. 20546. 


ii 



ACKNOWLEDGMENTS 

After a work is finished and one looks back to acknowledge the 
contributions, and express gratitude to the contributors, he finds it 
impossible to count everything and everybody. Thus, he is compelled 
to follow a conventional way, and mention only those who are closely 
related to that specific work, keeping in his heart the feelings toward 
the silent majority of contributors. 

The writer of this report wishes to express his deepest gratitude 
to his Supervisor, Dr. Ivan I. Mueller for the cooperation, understanding 
and support given fully and promptly always. 

Great indebtedness is owed by the writer to Dr. R. H. Rapp for 
the valuable suggestions and comments and the extra time he spent to 
make it possible for the writer to meet the deadlines. 

The writer feels grateful to the National Aeronautics and Space 
Administration for the financial support he was granted during his 
studies at The Ohio State University. 

Last, but not less, the writer wishes to express hearty thanks to 
the secretary, Mrs. Evelyn Rist, who, besides her regular duties, 
arranged for the writer to have this report typed, in spite of 
the restricted time limits. The writer is also thankful to Mrs. Irene 
Tesfai and Miss Barbara Beer for the excellent typing of this report. 

This report was submitted to the Graduate School of The Ohio 
State University for partial fulfillment of the requirements for the 
degree Master of Science. 


iii 



TABLE OF CONTENTS 


Page 

PREFACE. ii 

ACKNOWLEDGMENTS iii 

TABLE OF CONTENTS . . . . iv 

1. INTRODUCTION 1 

2. THE CASE OF TWO OBSERVING STATIONS 5 

2.1 Geometric Analysis. . 5 

2.2 Dual Adjustment 10 

2.21 The Duality Principle 10 

2. 22 Formulation of the Problem through the Duality 

Principle. . 16 

2.23 The Euclidean Geometry 18 

2.24 The Euclidean Metric 21 

2.25 The Usual Euclidean Metric and Its Dual. . . 29 

2.26 Statistical Analysis of the Adjustment Problem 47 

2.261 Statistical Model 47 

2.262 Comments on the Statistical Model 50 

2.263 Least- Squares Solution 54 

2. 3 Determination of the Cartesian Coordinates of the Stations 60 

3. THE CASE OF MORE THAN TWO STATIONS 60 

BIBLIOGRAPHY 64 


iv 



1. INTRODUCTION 

The science of geodesy has found in artificial satellites new fruitful means 
toward accomplishing its objectives. Thus the new branch of satellite geodesy 
was born, comprising the methods of treatment of the geodetic problems which 
employ the new celestial bodies. 

Satellite geodesy methods may be considered as falling in three categories: 

(a) geometric methods , employing mostly geometric principles, 

(b) dynamic methods , whose results rely mostly on dynamic principles, and 

(c) mixed methods , when no discernment can be drawn between the previous 
two. 

The above categories furnish the basis for the distinction between geometric 
satellite geodesy and dynamic satellite geodesy. 

Given a science, say geodesy, its objectives are more or less defined and 
what changes are the methods and the means scientists invent and employ to come 
closer to the solutions of the problems in question. Therefore, comparison of 
the results from different methods and means employed for the same objective 
is made very often. 

It seems reasonable to state that the more independent the employed methods 
and means are from each other, the more reliable are the inferences from the 
comparison of the corresponding results. 

Considering the case of geometric and dynamic satellite geodesy, there 
appear to be problems which may be treated with methods of either one of the 
two branches. One example is the determination of the equatorial radius of the 
earth. 

In this and analogous cases, one branch serves the other better, when the 
results rely as much as possible upon the principles of the branch within which 
they were obtained. Consequently geometric satellite geodesy problems should 



be treated in such a way that dynamic principles involvement is avoided as 
much as possible. 

It is illuminating to trace back to the fundamental principles upon which 
the geometric satellite geodesy is founded. 

The bodies of the earth and the satellite are considered in continuous 
motion in space. Their motions are not independent from each other, but never- 
theless we do not know ti c exact relationship between them because of dis- 
turbing factors of not well known effects. Geometric satellite geodesy wants 
to ignore any knowledge of the relative motion of the satellite with respect to 
the earth, and considers the satellite at certain discrete moments as a target 
point in space. Instead of the bodies of the earth and the satellite being in con- 
tinuous motion, now we have only the solid earth in continuous motion and a set 
of discrete points Q u Qa, . . . Q n corresponding to the satellite positions at 
times tj, t 2 , . i . t n . 

One should look at the problem as following. At time t t the satellite is 
at the point Q t while the solid earth has a unique relative position Ej with respect 
to Qj. This is an event which is given the designation TEi, Qj/tj], 

If there exists a means to snap, so to say, an event fEi.Qi/tjlin such a way that the 
relative position of Q r with respect to the solid earth may be recovered, then one 
says that the point Q t may be tied to the solid earth. Having the points Q t , Q 2 , . . . 

Q n tied to the solid earth we may consider the solid which comprises the solid earth 
and the finite set of points Q 2 , ... Q n . Let this solid be called extended field 
of geometric satellite geodesy. The field of the earth-tied geometric geodesy is 
defined on the set of points of the solid earth with operations the angle and distance 
measurements. The field of geometric satellite geodesy is an extension, so to say, 
of the field of the earth-tied geometric geodesy with the same operations, i. e. , the 
angle and distance measurements. Problems which could not be solved in the 
field of the earth-tied geometric geodesy are solved in the field of geometric 
satellite geodesy. This is a point where one is tempted to recall in mind the 
analogy with the algebraic fields. No matter how much apart in nature the two 
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cases are, one cannot deny the analogy of the above geodetic case with that of 
the fields of real and complex numbers. For example, the equation x +1 = 0 
does not have any solution in the field of real numbers, but it does have in 
the field of complex numbers, which is obtained from the former after the ad- 
junction of the element i = J- 1. If nothing else this analogy hints some beauty in 
formulation, deeper understanding and broad perception of the situation, gained 
when the concrete structures of the specific problem in question are identified 
with the corresponding abstract mathematical structures. 

A method of solution of a problem in geometric satellite geodesy might be 
considered, just for description purposes, as anticipating two things: (1) the 
recovery of the relative positions of the extension-points with respect to the 
solid earth, and (2) the solution of the problem itself. These two things are done 
together in a unified way; however, this remark is of importance for the formula- 
tion of the problem. 

The recovery of the extension point Q, of the event fEj, Qj/t^is achieved 
by performing simultaneous observations in a certain proper mode from a number 
of ground stations. Depending upon the observational mode and the problem to be 
solved, there exists a minimum for the number of participating stations at each 
event, and the number of events needed for a unique solution of the problem. 

Geometric satellite geodesy has at its disposition today three basically 
different observational modes: 

(a) the optical observations mode, 

(b) the range observations mode, and 

(c) the range difference observations mode. 

The problems of geometric satellite geodesy usually consist in determining the 
cartesian coordinates of a set of ground stations. It might happen some of them 
to be of known position. In this case the problem is easier than the problem 
where all the stations are of unkown position. Geometric satellite geodesy 
methods provide solutions to the general problem where all the stations are 
considered of unknown position. Geometrically speaking all the observational 
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modes are equivalent to each other with respect to that problem. 

In the optical observations mode each participating station P { at an 
event TEj Qj/tj] observes quantities that are geometrically equivalent to the 
directional cosines of the direction P t Qt with respect to a coordinate system 
which must be connected with the solid earth. It should be explained here that 
the coordinate system is considered connected with the solid earth either when 
it is tied to the solid earth, or when it is not tied but the relative motion is con- 
sidered known. Here and throughout this work the employed coordinate 
system is considered tied to the solid earth. 

The problem to be solved is the determination of the cartesian coordinates 
of the participating stations. The questions to be answered in the following 
sections are: 

(a) What information do the observations alone (without the introduction of 
additional data) provide about the configuration of the stations and what are actually 
the quantities to be adjusted ? 

'(b) What additional data are needed to be introduced in order to determine 
the coordinates of the stations? Are these additional elements necessary during 
the adjustment? 
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2. THE CASE OF TWO OBSERVING STATIONS 


2. 1 Geometric Analysis 

As it has been stressed earlier each method of geometric satellite 
geodesy must anticipate observations for the recovery of the extension points 
(satellite position points) on the one hand and the estimation of the unknowns of 
the pi'oblem on the other. 

Since two intersecting 
directions determine a point 
it is concluded that in the 
optical observations mode 
the minimum number of 
participating stations at 
each event is two 
(at least for the 
tying of the extension 
points Q t ). It will be 
examined presently 



Fig. 1 


Fig. 2 


whether that information at each event is enough to solve the posed problem. 

Let P r and P 3 Qj + l 

be two ground stations 
which observe simul- 
taneously at the instances 
rtj l - rtj , tg, • . . > tpl the 
satellite which at these 
instances is at the points 

r Qj ] = r Qi,Qs, . . . ,Q P ] 
respectively. 

At the instant tj the stations P t and P 2 observe the point Qj and obtain 
information equivalent to the directional cosines [ay. by, c u ] of the direction P x Qj 
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and ra 2j , b 2J , ] of the direction P 2 Q Jf with respect to some cartesian coordinate 
system defined above. 

TjOt the following notation stand for the adjacent discretion: 

X, , Yj , Zj Coordinates of Qj 

Xj, y t , z t Coordinates of ground station Pj 

ra 1Jf bjj, Cjj] = nij Directional cosines vector of the direction PjQj, where 

a, b, c correspond to the x, y, z-axes respectively. 
it j Plane defined by Qj and the two observed directions. 

There appear to be three distinct views, namely 

(a) the plane defined by the three points Pi, P2, Qj 

(b) the plane defined by the line PjP 2 and the point Qj 

(c) the plane defined by the intersecting directions PjQj and P^j which is 

the above designated by Try . v 

The question arises as to which one of these fits the actual situation in this 
problem. It is not difficult to single out the last one. As a matter of fact there 
are two observed directions toward the point Qj. Although the observations are 
made from P x and P 2 the station position is not specified on the respective directions. 
Thus one may not consider the points P lf P 2 or the line PiP 2 for the definition 
of the plane. There remains only n j and let it be called fundamental simplex of 
the optical observation mode. 

The equation of this fundamental simplex is that of a plane which passes from 
a point QjTXj, Yj.Zj] and is parallel to two vectors rHi = Tau.bu.Cul and m 2 = ra 2J , 
ta 2 j,c 2J ]. That is 


x-Xj 

•V-Yj 

z - Zj 



a u 

b lj 

C U 

= 0 , 

(1) 

a aj 

b aj 

C 2J 




or in normal form 

Qfjx + 8jy + y s z + 5j = 0, (2) 
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where 

b 1.1 c u 

. (3) 

b 2J . C 2J 
C U a lJ 

ft = ' . (4) 

c ag a a.i 

?u l>u 

Vj = . . , (5) 

a z b su 

Yj Zj 

a lj b u C 1J 

a aj b ai c sa 

Suppose that index j goes over the values 1, 2, . . . , p. Then there is a set of 
p planes it s (j = 1, 2, . , p). Notice that each of these planes is defined independently 
from the others. Any pair of these planes uniquely defines an undirected line. 

There exist C\ ~ p(p- l)/2 such lines. 

Consider any two from the planes , say j = j and j = k, i. e. , 

*r.i = + fty + Vjz + 5 3 = 0 (7) 

and 

TTj : OfeX + &y + Y,,z + 6k = 0. (8) 

Let tjic be the line of intersection of these two planes and a Jk , b 3k , c Jk its 
direction cosines. Now the following theorem of analytic geometry will be proved. 
Theorem The direction cosines of the line of intersection of two intersecting 
planes are proportional to the two-rowed minors of the coefficient matrix of the 
equations of the two planes, taken alternately with the plus and the minus signs. 

Proof. Consider any two points P t and P 3 of the line l Jk of intersection. Then 
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«jXl + 

ftyi + yj z i + ft 

= 0, 

O^Xi + 

Ayi + + 6, 

«,X 2 t 

ft y 2 vjz 2 • 5j 

= 0 , and 

0^X 2 + 

fty 2 + v k z 2 -< 5, 


»j(Xi - X 2 ) + 0 J (y 1 - y 2 ) + V 3 (z 1 - z 2 ) = 0 

and 

QtvX x - x 2 ) + ft(Yi " Ys) + yk(Z! - zg) = 0. 


But 


and then 


*l..--Xa Yi- Y2 = zi - z g 

a Jk ftk C jlc 

a j a jk + ftftk + VjCj k = 0 
a k a Jk + ftftk + VkCjk = 0 


A 


a Jk 


~ y i 

Ofk ft_ 


_ftk_ 


_7k 


By Cramer's rule the last system yields: 



ft 


a J 

-Vl 

^k 


• Cn and bn. = 


^k 

V J 

ft 

^Jk ouu u Jk 


ft 

A 

ft 


% 

ft 


a 1k _ ftk _ c Jk 


ft 

Vj 


ft 

ft 


ft 


ft 

Vic 


7k 

«k 


a k 

ft 



( 9 ) 

Q.E.D. 


It is of interest to find the value of 
Equation (9) yields: 
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1 


or 


l Jk 


1 ftv 1 C tk 


Jk 


ft Vj 

2 

-1 

V 

2 

-1 

ft 

2 

ft v k 


% «k 


Ofc ft 



ft 

yj 

2 

+ 


2 

+ 

ft 

ft 

Vk 


?k «k 


Ofc ft 


( 10 ) 


Therefore the direction cosines of i Jk are completely defined from the two planes, 
provided they intersect. 

As it was mentioned before there exist C| = p(p - l)/2 such intersection 
lines i Jk . The adjustment problem becomes obvious: What is the line cbest 
fitting the p(p- l)/2 lines L, k ? Statistically speaking, this is a regression 
problem where all the components of the position vector of each "point" are 
subject to error ( [22] pp. 186-194). 

The question springs up as to which quantity one should apply the least 
squares optimum criterion in such a regression problem? There exists a uniquely 
defined quantity between any two lines, namely their shortest distance which is 
along their common perpendicular; but this is a length, and there exists no length 
in the problem. 

The observations are angular quantities that means absolute invariants for the 
metric geometry and an adjustment of purely absolute invariants through relative 
ones, as the "length" is, would be considered a spoiled one. It seems reasonable 
to keep an adjustment as "pure" as possible. Therefore by principle the introduction 
of any "length" is excluded here. The writer drafts here the principle of duality, and 
prefers to go back enough into the fundamentals of geometry to found the adjust- 
ment to be employed for this problem. 

The following subsection is devoted to the application of that principle to the 
problem in question, i. e. , the best fitting line to a given set of lines. The result- 
ing adjustment drawn through the duality principle will be called dual adjustment . 
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2. 2 Dual Adjustment 


2.21 The Duality Principle 

A refresher introduction of the duality principle stressing the points that 
the problem in question needs more, is set forth now. Duality is not a theorem 
but concerns theorems. As any other principle in mathematics, duality principle 
belongs to the jurisdiction of "metamathematics. " For wider comprehension of 
the concept, it will be given from two instructively different standpoints: 

(1) that of synthetic geometry which studies figures without employing formulas, 
and (2) that of analytic geometry which establishes a correspondence between 
the elements of a figure (whatever they are considered) and the elements of a set 
called coordinates, and studies geometry through formulas between the coordi- 
nates corresponding to geometric relations. The intrinsic discussion will be 
within the latter point of view. 

Duality in Synthetic Geometry. 

Each geometric figure may be considered as an assemblage of elements, i. e. , 
"units" from which the figure may be constructed in some way. For example a 
plane-ellipse might be considered as the aggregate of points with which one of 
them moving along the ellipse succesively coincides. But equally well that plane- 
ellipse might be considered the envelope of the coplanar tangents at each of its 
points. Thus the same figure has been generated by using two different elements, 
i. e. , the point and the unlimited line in the plane of the ellipse. This is the case 
for any plane figure. Analogously in space each figure may be considered as 
generated either from a moving point or from a moving plane. It is this double way of 
viewing a figure being generated on which the principle of geometric duality is based. 

It was Poncelet, who first enunciated the principle of duality, which in the 
case of the plane is: 

Any theorem about properties of position of plane figures is 
accompanied by the so-called dual one, which has the words "point" 
and "line" mutually interchanged. For example the dual statement of 
"three points in a plane define a triangle", is "three lines in a plane 
define a triangle. " 

The duality principle in space is: 
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Any theorem about properties of position of figures in space is 
accompanied by its dual one which has the words "point" and "plane" 
mutually interchanged, while the word "line" in place. For example 
the dual statement of "three not collinear points define a plane, " 
is "three not coaxial planes define a point. " 

In the above statements of duality there is no mention of the so-called mass 
relations, i.e. , those which include distance angle, area, etc. In faet, the 
duality principle was introduced within the frame of projective geometry at a 
time when it was unknown how to deduce metric geometry from projective 
geometry. * The Cayley's principle that metric and affine geometry can be con- 
sidered special cases of projective geometry made possible the extension of 
the duality principle over the mass- relations of metric geometry. This is the 
key for the solution of the problem in question. The analytic character of this 
problem naturally turns the discussion into the field of analytic geometry. 
However, the fact that the whole deal is in the three-dimensional geometric space, 
that means within visualization, helps somehow against unrealistic analytical 
representations. 

While in synthetic geometry duality is based on the double way of viewing 
the geometrical figure being generated, in analytic geometry duality is based on 
the double interpretation of the formula that stands for a figure or its analytical 
representation with respect to a coordinate system introduced beforehand. 

Duality in analytic geometry is extended in a unified way to spaces of any 
dimension. However, this discussion need not go beyond three-dimensional 
space. Some representative examples will now illustrate the principle. 

For the duality in plane consider the equation of a straight line in the plane with 
respect to a cartesian coordinate system, (orthogonal or oblique) i.e. , 

ax + by + c = 0 (11) 


* The English geometer Cayley in 1859 presented a way to deduce affine and metric 
geometry from projective geometry. 
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In homogeneous coordinates x h x 2 , x 3 related to x and y through the relations 
x = and y = ^ , (x 3 i 0) 

x 3 x 3 

equations (11) may be written as 

ax! + bx 2 + cx 3 = 0. (ll/ 

The last relation is a homogeneous linear equation with respect to x L , x 2 , x 3 and 
the same holds with respect to a,b,c. One could say that the relation is symmetric 
with respect to the triples {x lf x 2 , x 3 | and<{a, b, c| . 

Now to what extent are the a, b, c, conversely determined by the 
straight line ? If one compares (11)* with, say, 

ax t + bx 2 + cx 3 = 0 (ll) 77 


under the assumption that they represent the same straight line one has, considering 
the linear system (11)* and (11)", that * 


which implies 




Xc*. 


( 12 ) 


Consequently, the quantities a,b,c are determined by the straight line up to a common 
constant of proportionality. 

Similarly each point determines x u x 2 , x 3 up to a common factor of pro- 
portionality. Indeed, if (x/, x 2 *, x 3 *j represents the same point with {x lt x 2 , x 3 | 
whose non-homogeneous coordinates are x and y then 


/ ■ • / 



* Since r a, b, c] and Ta , b , c ] are two non- vanishing linearly dependent vectors, 
each is a multiple of the other. 
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or equivalently if x 3 = 71 x 3 , then x t = x • X 3 = xfix 3 = u(xx 3 ) = fjiXy and 
x 2 = y Xg - y * MX3 = |i(yx 3 ') •= uxb, i.e.. 


M x 1 , x 2 


M x s < 


M x 3 


(13) 


The important fact here is that one can associate coordinates to the straight 
line of the plane in the same way as with the points of the plane. Consequently 


each triple 77 , may be interpreted either as a point or as a straight line 
of the reference plane. Now, all the theorems that state anything about the 
points of the reference plane, about the aggregate of such points, and about the 


£,t?,C , 


the 


relations among them are nothing but statements about triples 
aggregate of such triples and algebraic relations among them. But the very same 
triples may be interpreted either as points or lines of the reference plane. 

Two statements, one being expressed in point coordinates and dealing with 
points, the other in straight line coordinates and dealing with straight lines 
correspond exactly to each other and are called dual statements . Analogously 
"point" and "line" are called dual elements in plane. As one realizes, dual 
statements are different interpretations of one and the same algebraic result. 

In order to have a geometrical visualization of this example, that synthetic 
geometry would provide, one could proceed this way: Keep the straight line 

3. C 

constant in the plane (i.e. , — = a. = constant and — = j3 = constant) and let the 


point | x = -^v Y 
points of the line. By keeping the point < x = 


— 2 ; varying. Then equation (11)* represents the range of 

Xg J 

I x, Xo'l • 

y = — * r constant in the plane, 

X»a J 


x 3 


and having the straight line — - cl, — = )3 (-varying the very same equation (11) 

^ 3. E J 


represents the flat pencil of lines through the point x 


2b 

x 3 


y = 


X, J 


For the duality in three-dimensional space consider the equation of a plane, 


say, 


ax + #y + yz + 6 


0, 


(14) 
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or in homogeneous point coordinates x t , x 2 , x 3 , x 4 irelated to the non-homogeneous 
x, y , 'h through the relations 



y 


-*2 

9 

X4 


^ , (X4 / 0) 

A4 / 


ax, + /3x 2 + vx 3 + 5 X 4 = .0. (14)' 


Relation (14)' is a linear homogeneous equation with respect to x^ x 2 , x 3 , 
x 4 and the same holds with respect to a, $, y, 6. 

In the following; the same steps which were followed for the equation of a straight 
line in a plane above are repeated. In fact, both of these cases are special cases of a 
general theory about hyperplane coordinates of a general projective space. Under 
the general theory, the commonly conceived planes are the "hyperplanes" of the 
three-dimensional projective space and the straight lines are the "hyperplanes" 
of the two-dimensional projective space. Nevertheless the same steps are re- 
peated, for this is the space (i.e. , the three-dimensional projective space) to which 
the problem in question belongs. 

First, the question: to what extent are a., jS, v, 6, conversely deter- 

mined by the plane? To answer this question equation (14)' is compared with, say, 

O'Xj + jS'x 2 + v'x 3 + 6x4 = 0 (14) " 


under the assumption that they represent the same plane. Then 

fa (3 y 5~| 


rank 


= 1; 


[os 0', v', 5'J 

since the intersection of (14)'and (14)" is just the plane itself, that means of 
two dimensions, and consequently these planes are linearly dependent. Therefore, 
there exists X/ 0 such that 


ot = \ct, j3 = .AjS', v = Xv', 6 = X5'. (15) 


An alternative procedure is to write 
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(orx! + 0x 2 + vxg + 6 x 4 ) - X(a"x 1 + 0x 2 + y x 3 + 6 x 4 ) = 0 
or 

(a- X.Of'jXi + (0- X0')x 2 +'■ (v- Xy')x 3 + ( 6 - X 6 ')X 4 = 0, (16) 

for the linear dependence of the planes; then from the fact that (16) must be 
identically satisfied, (15) follows. Consequently, the a, 0, y, 6 are determined 
by the plane up to a common constant of proportionality. 

Similarly each point with non-homogeneous coordinates x, y, z determines 
x i> x 2 , x 3 , x 4 up to a common constant of proportionality through the relations 



Indeed if x/, x 2 / , x 3 , x 4 * represents the same point with x,, x 2 , x 3 , then 


■ • *1 = JL v - *z = 2% z = '2Ea. = 25a! 

x ■ — f , y — * , h — , ' , 

X4 X 4 X4 X4 X4 X4 

or, by putting X 4 - /ix 4 it is obtained 

x x - fj, x x 2 =-'jLi.Xa # , x 3 = 11X3', X4 = n*. 4'. (17) 


Thus, one can associate coordinates to the plane in the same way as with 
the points in space. * Consequently each quadruple Jjf, rj, C, rjmay be inter- 
preted either as a point or as a plane. This is the point where duality principle 
in space springs up. For, all the theorems that state anything about points in 
space, about aggregate of such points, and about relations among them are nothing 
but statements about homogeneous quadruples j^, T), C, tJ, aggregate of such 
quadruples, and algebraic relations among them. But the very same quadruples 
may be interpreted either as points or planes. Two statements, one being expressed 
in point coordinates and dealing with points, the other in plane coordinates and 
dealing with planes are called dual statements . 


* Here by space it is meant the three-dimensional projective space. The adjective 
projective has been dropped to avoid confusion for some readers. 
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Again one notices that, dual statements are different interpretations of 
one and the same algebraic result. The above example exhibited that "point" 
and "plane" are dual notions in space. No mention was made about the third 
fundamental element* of space, the straight line. But the dual statements 

"two not -coincident points define a line" 

"two not-coincident planes define a line" 

show clearly that "straight line" is a self-dual notion in space. 

After the above general introduction of the duality principle, one has all 
that is necessary to formulate the problem in a way that turns out to be very 
fruitful. The following paragraph is devoted to this formulation and constitutes 
the key for the adjustment problem. 

2.22 Formulation of the Problem through the Duality Principle 

The geometric analysis of the optical observations mode from two 
observing stations (section 2. 1) showed that the adjustment problem 
amounts to that of determining the best fitting straight line of a given set of 
straight lines in space. Each straight line of the set to be fitted is the intersection 
of a pair of planes from a finite set of planes. Therefore the adjustment problem 
may be stated: What is the best line (in some respect) through a set of planes ? 
Instead of attacking this problem directly, it is easier to consider its dual one, 
which is a two-dimensional problem. The formulation of the dual problem is obtained 
by interchanging the words "point" and "plane" and leaving the word "line" in place. 
Thus the enunciation of the dual problem is: What is the best line (in some respect) 
through a set of points? This is a problem in two dimensions (see section 
2. 21), and its analysis occupies a first place in regression analysis. The 
optimum criterion for the ,r best" will be that of least squares. What is needed 
from this problem is the procedure of solution and not the solution itself. 

*The characterization fundamental element for "point", "plane" and "straight line", 
does not have any geometric justification. It is given on intuitive grounds only; 
for there are other "elements" with equal geometrical justification. 
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Therefore, the solution procedure of the two-dimensional dual problem will be 
sketched first, and then it will be translated through duality for the problem in 
question. The problem of fitting a straight line to a set of points in a plane 
depends upon how the coordinates of the points are treated in the adjustment, 
from a statistical point of view. Thus there exist two cases: 

(i) Only one of the coordinates of each point (either x or y) may be subject to 


error, and 

(ii) both of the coordinates may be subject to error. 

In the first case, one obtains either the regression line of y's and x's (when 
only y's may be subject to error) or the regression line of x's on y’s (when only 
x's may be subject to error). 

In the second case, under some assumptions, there exists a unique regression 
line. The case of interest here is the second one. As it is very well known, the 
regression line in this case results from the minimum of the sum of the squares of 
the distances of the points from the line to be fitted, provided the two coordinates 
(x and y) are given equal weights*. If the coordinates x and y are given equal 
weights, then the line of the observed and adjusted point is perpendicular to the 
fitting line. It is supposed that this is the case here. 



*If the coordinates x and y of a point P are not given equal weights, then the line 
connecting the observed and adjusted points is not perpendicular to the fitting line. 
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LetPjfx^y^be a set of points in aplane and (e) the least-squares line through them 
(Fig. 3). If P 0l is the adjusted point of P^ then (PiP ol ) r (e), and the least-squares 

criterion is S df = min. 

1 

Nowlet('7r 1 )be asetofplanes in space, and (£) the least-squares line through them 
(Fig. 4). Then duality principle establishes the following correspondence between the 
elements of Figs. 3 and 4. 

(O — (0 
Pol — ri ol 
(PoiPi) — » m 
(0 X (PolPl) — - (t) x(4). 

Let D(.ir ol , TTj) be a symbol standing for a concept of "distance" between the planes 7r 0l and 
rij. Then the least-squares criterion £ df = min. forpoints is translated into 
£ lfi(ir oi ,TT) = min. forplanes. Thequestion now arises as towhatis the "distance" 
between two planes. One might think of the angle between the planes. At least intuitively 
angle is a suitable measure here, but there are many ways to measure a quantity, which 
are equivalentin some respect. Therefore.it is not hopeless to search for a suitable 
convenient measure of the "distance"between two planes. This is the target the 
following subsections aim at. 

2.23 The Euclidean Geometry 

The problem in question falls within the content of the three-dimensional 
Euclidean geometry. Although the Euclidean space is intuitively appealing, there exist 
some questions which are of importance not only theoretically, but also from the 
point of view of the applications. 

It is not the purpose of this work to go over the foundation of the Euclidean 
geometry. However, trying to define an appropriate metric for the adjustment 
problem as set up above, it is inevitable to face the following question: What is con- 
sidered a legitimate metric in Euclidean geometry and how is one to determine the 
totality of these metrics? This question does touch the foundation of Euclidean geometry. 
A rigorous presentation entails the parallel exposition of non-Eucli dean geometries in the 
content of which Euclidean geometry is a special case. In order to reach the problem's 
goal without the employment of non-Euclidean geometries, for that would take the 
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discussion too far, a conventional way will be followed, where rigorousness will be 
guaranteed by recalling theorems of the general theory. 

In any system of geometry the existence of some fundamental elements, which 
constitute the objects of the operations of the geometry, is assumed. Different 
kinds of geometry arise considering different fundamental elements or by defining 
them in different ways. From this point of view in Euclidean geometry the existence 
of two objects is assumed: 

(a) the generating element (either point, or straight line, or plane, or some 
other configuration of equal justification), and 

(b) the distance between the generating elements. 

Of course these are only the objects. In order to construct a geometry a certain set 
of axioms is needed upon which the structure of the geometry will be based. If now 
someone takes the point as the generating element and the distance between two 
points in the commonly known sense and adopts the set of axioms of Euclid, the 
point-Euclidean geometry is obtained. If instead of the point, one takes the straight 
line or the plane as the generating element, one obtains equally justified geometries: 
the line-Euclidean geometry or the plane -Euclidean geometry, respectively. 

Geometries different from the Euclidean may be obtained either when the 
distance is defined in a different way, but equally justified (in the sense that a 
self-consistent geometric structure* may be built by using that), or when there 
exist different, equally justified sets of axioms. The latter case is that of non- 
Euclidean geometries, which arose by disputing Euclid's fifth axiom, the so-called 
parallel axiom **. 

From the above point of view, geometry seems to possess no organic 
systematization. This, however, is not the case. Analytic representations in 
geometry appeared to be very fruitful and enabled geometers to develop geometry 
in a systematic way as a unified structure. 

*One may recall here as an example the Minkowskian geometry versus Euclidean 
(see [11], p. 133). 

**Some geometers distinguish between axiom and postulate, the former being a 
simple fact of logic while the latter is about space perception. 
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It was Klein's ingenious idea to introduce the notion of "motion" of figures 
in geometry which has provided the best systematization of geometry so far. 

Klein's view of geometry is simple, elegant and fruitful, and it will be adopted 
here. 

In analytic geometry one represents the various geometric figures in terms 
of coordinates. For example, one thinks of the totality of points in the usual space 
as represented by the totality of triples x,y,z, i.e. , by their coordinates with 
respect to some coordinate system. To every transformation of points in space 
there corresponds a certain transformation of their coordinates. The converse is 
not true. That is, not every analytic relation between coordinates corresponds to 
a geometric transformation. In other words, each system in geometry (or each 
"geometry" as it is usually said) possesses a well-defined group of motions, which 
can be represented analytically (in terms of coordinates). One defines then the 
geometry as the theory of the properties of the configurations which remain un- 
changed under those motions. In this light. Euclidean geometry is defined as the 
geometry which studies those relations between the coordinates, which remain 
unchanged by the linear substitutions of coordinates corresponding to the following 
motions: 

parallel displacements, 

rotations about the origin, 

reflections about the origin, and 

similarity transformations with the origin as center. 

These substitutions are linear and constitute a group of transformations, the so- 
called principal group* of transformations . For a figure to remain unchanged 
under these "motions" means that it is "moved" as a solid in space; that in turn 
implies distance and angle preservation. 

In the principal group, the motions involve six parameters, to which one must 
add one more parameter for the change in unit length (scale), so that altogether 

*This set of transformations constitute a group, for (i) the product of any two of 
them belongs to that set, (ii) it comprises the identity transformation and each 
transformation has its inverse one. 
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the principal group needs seven parameters to be defined completely. 

Let this group be designated by G 7 . 

Geodesists are very familiar with this group. In each geometric adjust- 
ment in geodesy, the conditions and the constraints involved arc dictated by the 
principal group as it is considered in connection with the geometrical configura- 
tion and the coordinate system of that adjustment. 

2.24 The Euclidean Metric 

Consider the three-dimensional geometric space, * as the set either of its 
0 points or its straight lines, or its planes. Let E|, E 3 , be the designations 
'for these three sets respectively. These are the three sets of interest of the 
problem in question. As it was shown earlier (while discussing the duality 
principle) one may introduce a coordinate system, and establish a correspondence 
between the elements of the sets E£, E 3 , and the coordinates. Specifically 
there was established the correspondence of 

points with the quadruples r Xl> x s , x 3 , x*], 
planes with the quadruples fa, 8 , y, 5 ], and 
lines with the pairs of points or planes. 

It was said also, that each point determines the quadruple fa^Xs, x 3 ,x 4 ] up to a 
common constant of proportionality, and similarly each plane determines the 
quadruple fa, 8 , V, 61 up to a common constant of proportionality. That is to 
designate 

r Xl , x 2 ,X3,x 4 ] ; S rxx t , Xx 2 , XX3, AX4! and 

fa, 8, V, 6] - r\o;, >8, Xv, X5]. 

Remark. If x 4 / 0, one may write 

X = ^ y = *2 = 2*2 z = *3 = 2*2 . 

X 4 Xx 4 X4 Xx 4 x 4 • Xx 4 

♦This is the empirically conceived three-dimensional space. 
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But what does the condition x 4 ^ 0 really mean? One can easily answer the 
question by taking x 4 -»0. Then x -»°°, y ->°°and z So the condition X 4 ^ 0 
excludes the points at infinity. Since the configuration of the ground stations and 
the satellite positions is considered of finite extension, one can exclude the 
points at infinity without any consequences for the problem. When x 4 ^ 0 the 
convenient value x 4 = 1 is usually chosen, and the quadruples l’x 1 ,x 2 , x 3 , 1] are 
called affine coordinates. One then observes that x = x t , y = x 2 , z = x 3 . Never- 
theless, x* will be kept in the formulas for homogeneity and symmetry purposes, 
and it will be always assumed in the sequel that X 4 ^ 0. One might ask, whether 
it is necessary to exclude the points at infinity or not. The answer is not 
supposed to be given here, but anyway it is in the affirmative, for distance is 
not defined at infinity. 

Going back after this digression, one observes that, points and planes 
have the same analytical representation. This is very important, because it 
allows for simultaneous analytical treatment of both the two sets E 3 and 
without having to distinguish between '’point’’ and ’’plane’’ during the analysis. 

For a simultaneous analytical treatment of the sets E 3 and e£, a general set E 
is introduced, which designates the set of any Euclidean space to be defined 
below. 

Having the non-empty set E, the notion of distance between the elements 
of the set naturally suggests itself. The non-empty set E equipped with a suitable 
concept of distance, (provided there exists a suitable one), is called Euclidean 
space. But what does one mean by distance, and what does this notion render? 

In the set E, a certain type of convergence is defined for the elements of the set. 
For example, in the sets E 3 and E^ the following types of convergence are defined 
(see Figs. 5 and 6). Let rp 0 , P lt . . . , P 1( . . . ] be a sequence of points in E 3 . 

Then one says that this sequence converges to the point P, if the representations 
of the points P^x^ y lf z t \ i - 0, 1, . . . , i, . . . , converge to the representation of 
the point prx,y,z], i.e., 
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r Xo,X v , • » X 1. • • ■ 1 ,-> x 

ry 0 ,yi, . . . ,yi, ... y 

r Z 0 , Z^, . . . , Zj, . . . ] z 



Similarly, if n 0 , v lt . . . , n lt ... is a sequence of planes in Ej, then one says that it 
converges to the plane tt, if the representations of the planes rr^at, j8t, >i, 6il, i = 0, 1, 
i, . . . , converge to the representation of the plane rr^a, 0, v, 6], i. e. , 


Tofo, oil, .... oil, a 

r Bo, jSj, . . . , J3j, . . . 1 -» £ 

rv 0 , y lt . . . , v lt . . . 1 ■-* v 

ra 0 » ...1 

Note: The above Figs. 5 and 6 illustrate two special types of convergence in 
E3 and respectively, which are the only ones of interest in the problem in 
question. Fig. 5 illustrates the convergence of the sequence of the points 
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TPo, Pi, . . . , P lt ... 1 to the point P along the straight line of the points P 0 
and P. In other words it illustrates the convergence of the sequence of the points 
which is generated by P 0 , while it is approaching P along the straight line 
(P 0 , P). Fig. 6 illustrates the convergence of the sequence of the planes 
p!r 0 , tTl . . . , tTl •••"', which is generated by rr 0 while it is "approaching" the 
plane v , by turning around the line (-t) of the intersection of the planes tt 0 , and it. 
Now if there exists a function d of the coordinates of the elements of the set E 
which converges in an equivalent way to the convergence already defined in E, 
then d is called distance or metric . The set E along with the metric d, constitute 
a metric space designated by (E,d), which is the so called Euclidean space. 

The question faced now is: does a suitable metric d exist in the set E such 
that convergence in the metric sense is equivalent to convergence in the above 
defined sense in the set E ? It is reminded that E stands for either of the sets 
and Eg and since they have the same analytical representation their metrics 
(as functions of coordinates) will be of the same analytical expression interpreted 
in a dual way. 

At this point a result of the general theory about the foundation of Euclidean 
geometry is recalled, in accordance with what was said at the beginning, about the 
discussion of this topic. Thus the general foundation theory of Euclidean geometry 
proves: The general expression of a Euclidean metric d(P 0 , P) in three-dimension 
Euclidean space is the non-negative real function 


d(P 0 , P) 


l,k=l 


A x i 


*k _Xk 

X4 x^ 

. 

X 4 x 4 


(x£. x 4 ) 2 J 


(18) 


wnere [x°, Xg, x§, x 4 ]and [x L , Xg, x 3 , x 4 ]are tne analytic representation of tne 
elements P 0 and P. One may think of P 0 and P either as points, or planes, or 
any other equally justifiable configuration, with the same analytical representation. 
In non-homogeneous coordinates, expression (18) has the form 
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d(P 0 ,P) 


3 


I 


gik ( x ° - X i){^ - Xk) 


1 1 k=l 


(18)' 


which is obtained by choosing for the non-vanishing the value of one. 

Although (18) 'is simpler than (18), the latter is more useful for the problem in 
question. Before specifying the expression (18) for the cases of interest here, 
(points and planes), a few general remarks about the Euclidean metrics are in 
order: 

(a) The square of the Euclidean metric d(P 0 , P) is a positive definite quadratic 
form. Expression (18) may be written in matrix notation as follows. 


d 2 (P 0 , P) 




4 

x i 


4 

Xk 

4 

X 4 

• 

X 4 

X 4 


- V’ 


Xp.p G X, 


p„p 


(19) 


1>=1 


(x° • x 4 ) s 


where 


G = g* 


gn gis gia 

gsn gas g23 

gbi ga» g33 


G = G t 


( 20 ) 


and 


X, 


T _ 
P n P 


x ° X 1 


X 2 

Xg 


4 X.3 

4 x 4 


X° 

X4 


x° x 4 


X4 


X 4 ■ X, 


x° 

4 a 4 


X4 


? 0. 


( 21 ) 


In the case of non-homogeneous coordinates, (21) has the form 

x; p - r ( X °1 - ), (xg - x s ), (X§ - Xg )] 

= [ (x° - X), (y° - y), (z° - z)] . (21)' 

One says that the quadratic form d 8 (P 0 , P) = X T GX is positive definite, when 
X T GX > 0 for every X ^ 0. The same can be expressed in terms of the 
matrix G. Thus X’GX is positive definite when G is positive definite, meaning that 
it has positive eigenvalues. Since G is a real symmetric matrix, if X s , Xq 
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are its eigenvalues with p 1( ps, Pa the corresponding normalized* eigenvectors, then 

P t GP =■ A =. diag ( Xi, X a , A 3 ); (22) 

where P = p 3 , p 3 ]is the matrix whose columns are the normalized eigen- 

vectors of G, and A is a diagonal matrix whose diagonal elements are the 
eigenvalues of G. It is a well-known result that P is an orthogonal matrix, and 
P T P = PP T = I. Also P T - P\ 


(b) In the previous subsection it was said that the motions of Euclidean geometry 
preserve the distances and angles. The converse is also true. Having now both 
Of these measures (distances and angles) expressed analytically under one expression 
d(P 0 , P), one can state that the motions of Euclidean geometry are those which 
preserve the value of d(P 0 , P) for any pair of elements (P 0 , P), provided none of 
P 0 and P is at infinity. Let S be a matrix of a Euclidean motion. It is interesting 
to find out what kind of matrix S is. Consider the "distance" d(P 0 , P). Under the 
motion S, d(P 0 , P) is by definition an invariant quantity. If Y PoP is the new representa- 
tion vector, into which X PqP is transformed under S, one has 


Then 


Y PoP = SX PoP . ' (23) 

d(P 0 ,P) - x; oP G X PoP = (S 1 Y PoP ) T G (S 1 Y PoP ) 

= y; oP [(SYG(S^)] Yp oP . 

- Y ? 0 P MY P 0 P 


where 


M = (SY G(SV 




(c) One point of importance here is the following. In the case of a general metric 
space, analysis defines a metric d on a set E* as a real function of ordered pairs 
of elements of E which satisfies the following three conditions: 

(i) d{e lt e 2 ) £ 0, and d(e 1 , e 2 ) = 0 if e t = e a ; 

(ii) d(e lt e g ) = d(e 3 , e^, (symmetry); 

(iii) d(e t , eg) ^ d(e 1 , e$) + d(e 3 , eg), (the triangle inequality). 

The above given Euclidean metric obviously satisfies the conditions of this 
definition. Although on one and the same set E, one may define many metrics in 
the above sense-, each one renders different metrization of E. As it was said, 
the introduction of a metric follows the definition of a certain type of convergence 
for the sequences of the elements of the set. It is the type of convergence defined 
on the set which constitutes the criterion for the suitability of a metric. 

The Euclidean metric given by (18) is one of the so-called projective metrics 
The metrics of all the metric geometries derived from projective geometry are 
characterized as projective metrics. Thus, the use in some problem of a metric 
whose form is not that of (18) means that the problem is not treated within the 
Euclidean geometry goemetry. For example, the metric 

A x t 

di (P 0 » P) = max x° x 4 , (x°, x* i 0), 

1 = 1 , 3 , 3 

(< X *) 3 

or for x° = x 4 = 1 

di(P 0 , P) = max \x° - x t | (25) 

i=i, a a 

does satisfy the conditions of the above definition of a general metric, but it is not 
a projective metric. The same with the metric 

3 

d 8 (Po.P) - y Xl Xl (x°, x 4 ^0) 

Lj „0 X 4 X 4 

1=1 x 4 -^4 

' • ' • . \ 

*Here E stands for an abstract set. 
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or for x° = x 4 = 1 


d3(P 0 .P) = £|xg - x t | . 


(26) 


1=1 


(d) In this note the corresponding general expression for a Euclidean metric in 
two dimensions will be given. This can be done immediately by omitting the 
coordinate x* and considering the remaining coordinates x lt Xg, X 3 as homogeneous 
in two dimensions, i. e. , 


Then 


X, 

x = 

xa 

y =•-*, 

x 3 ¥■ 

0. 

X 3 

X 3 

3 

X° Xl 


d(P 0 , P) = 

I Su* . 

xg X. 

xg Xa 


1 ,k=l 

(4* 

X3) 3 


or in matrix notation 

rT 

'p 0 p ^p 0 p 


d 2 (P 0 ,P) = Xj oP G a Kp p, 


where 


G2 - (Slk) 


gu Bia 
S21 SSB 


and 


X*p = 


xl Xi 
x£ xg 


Xg Xg 
X 3 Xg 


xg • x, xg • Xg 


Gg - Gg, 


, X 3 / 0. 


In non-homogeneous coordinates, (27) and (28) have the forms 


* 

d(P 0 , P) = .[£ gik (x? - x t )(x£ - x*)] 2 


l,k=l 


and 


XffeP = [(X? - Xi), (4 - Xg)] 

= [(X° - x), (y° - y)] 


(27) 


(27) 


(28) 


(27)' 


(28) 
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2.25 The Usual Euclidean Metric and Its Dual 

The general expression (18) of a Euclidean metric is recalled here 


d(P 0 , P) 



1, k=l 


4 x i 


4 

x? x 4 


4 *4 


(4 * X4) 9 


1 

2 


, (Xl, X4 / 0) 


( 29 ) 


where [x° lf Xg, xg, xg]and [x lt x s , X 3 , x 4 ] are the analytic representations of the 
elements P 0 and P. At this subsection the elements P 0 and P will be interpreted 
either as points or planes. The matrix formulation of the above expression is 
recalled also from (19), (20), and (21). 

d 8 (P 0 ,P) - Xlj, G Xp o p, (29)' 


G = 


gll gl2 gl3 
gsi gas gaa , 
g3i gas gaa 


G 




4 xi 


Xg X 2 


xg 

Xa 

4 X* 

* 

4 *4 

9 

4 

*4 


4 • *4 x 4 • X* 34 • x 4 


, (4, x 4 y o) 


(30) 


While discussing the principal of duality in the three-dimension projective 
space, it was shown that "point" and "plane" are dual notions. It is true that the 
duality principle does not hold in Euclidean geometry in general. Thus whenever 
duality principle is drafted in Euclidean geometry, its validity must be demonstrated 
sufficiently. That is why the preceding analysis has been placed before this point. 

If in the general expression of a Euclidean metric one takes the elements P 0 
and P as points, then d(P 0 , P) renders a means of "distance"* between tne points P 0 
and P. If one takes P 0 and P as planes, then d(P 0 , P) renders a means of "distance" 
between the two planes. In the following, the first case will be considered first, 
for it is simpler and more visual than the second, in order to find out how the 
usual Euclidean metric results from the general expression (29). Then the second 


*The quotation marks are justified by the fact that the word "distance" stands for 
the usual Euclidean metric, which, as it is proved in the sequel, is a special case 
of the above. 
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case will be deployed on the first, using the duality principle as a guide, to find 
out the geometrical meaning of the dual metric to the usual Euclidean metric. 

Case 1 . Point-Euclidean space (E|, d). 

In this case the elements of the space are the points* of Eg, which the human 
mind seems to visualize better than any other configuration. The usual Euclidean 
metric, i. e., the usual expression for the distance between two points, will be 
connected with the general expression of a Euclidean metric. 

Let P 0 and P be tv o points with Cartesian coordinates [x^, x%, x 3 ] and 
[x lf x 3 , x 3 ] respectively. Then the usual Euclidean metric has the form 

d(P 0 , P) ■= [(x° - Xi) 2 + (x|- x s ) 8 + (xg - Xa) 8 ] 2 . (31) 


Introducing homogeneous coordinates, the last expression may be written 



X° X! 

s 

x a x 2 

2 

X3 X3 

d(P 0 ,P) = 

x£ X* 

+ 

X° X4 

+ 

4 X4 


(4**) 

3 • 

(X4 X4) 

a 

(x^x^ 2 


(31) 


or 



“ - 



. - 
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0 
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d(P 0 ,P) = 

Xl oP 

0 

1 

0 

Xp 0 p 



0 

0 

1 



where 



i 

~1 2 
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A x 3 


X3 x 3 

Xp o P = 


4 x 4 

* 

4 x 4 


4 x 4 


X4 x 4 


x 4 x 4 


x 4 X4 


Thus the usual Euclidean metric results from the general expression of a Euclidean 
metric for G = I 3 , i, e., when G is the identity matrix. Recalling the relation (24) 
of the remark (b) which followed the general expression (18), we observe that the 
matrix S of a motion in Euclidean geometry satisfies the relation 

S S T = S T S =' la (33) 


* It is recalled that when an element is considered as the generating element of 
the space, it is supposed to be undefined for the geometry to be constructed 
with that. 
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That means that S is orthogonal. Thus the motions in Euclidean geometry with 
the usual Euclidean metric are analytically represented by orthogonal matrices. 
Analogously, in the two-dimensional Euclidean space, the usual Euclidean metric 
results from the general expression ( 27 )' for 


G2 - 


1 

0 



Case 2 . 

In this case the elements of the space are the planes of E 3 . Although they 
have the same analytical representation with points, it is much more difficult to 
visualize configurations of planes than of points. Therefore, duality contributes 
substantially to problems with configurations of planes like the present one. The 
dual metric of the usual Euclidean metric will be formed and its geometrical 
meaning will be given. 

Letff 0 andff be two planes with analytical representations To^, j3 0 ,y 0 , 6 0 ] and 
[a, j3, y, 6] respectively. Their equations are 

a 0 X! + &x 2 + y 0 x 3 + 60X4 = 0, 

and (33) 

ax t + /3 x 2 + yxa + 6 X 4 = 0. 

Since the points at infinity are not considered in the problem (i. e., x 4 4 0), one 
could write these equations as follows: 

%X! + 0 o x 3 + y 0 x 3 + 6 0 = 0, 

and (33/ 

OX; + f3x s + yx 3 + 6 = 0. 

They resulted from the previous ones for x 4 = 1. By taking now the elements P 0 and 
P in the general expression (29)* of a Euclidean metric, as the planes ir 0 andff 
one has 

d(nr 0 ,7T) = [(« > T G(0 (34) 

TT 0 TT ff 0 ff 
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where 


(ft ) T = 
rr 0 ir' 


«0 

a 


So 0 


y 0 v 

$0 

6 

j 

6 0 6 

> 

6 0 6 


6 0 6 


6 0 6 


6 0 6 


(60, 6V 0) 


(35) 


and G remains exactly the same, for it is independent of the given interpretation 
to the representation vector Xp oP . Now for 
1 0 0| 

= I 3 , 


G = 


0 10 
0 0 1 


relation (34) yields 

d(7To,ff ) = un- ) T (ft„ )f. 

TTofTT TToiTT 


(36) 


This metric is a Euclidean metric, and, consequently, a projective one. It 
expresses the "distance" between tt 0 andtr in the same way the metric d(P 0 , P) in 
(32) expresses the distance between the points P 0 and P. But what is the geometrical 
meaning of saying "the plane TT "? Any comment on this question is postponed 
until later. Now the geometrical meaning of d(ir 0 ,TT ) is in order. While trying to 
interpret d(ff 0 ,fr ) geometrically, the above question will manifest itself in a 
natural way. 

Let 

Po = [o«f + 

and 

p = [a 8 + fi 3 + y 8 ]^. 

) T may be written as follows. 


Then the expression (35) for .(ft 
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where both terms of each fraction have been multiplied by l/(p 0 p). Since d (no,n ) 
is an invariant under Euclidean motions(see section 2. 23), a parallel displacement 
is applied on the configuration of the two planes n 0 and it such that 



Po P 

This transformation, though applicable, is not very clear from the present 
point of view. Therefore, another proof, purely analytical, will be given. The 
proof to follow, though rigorous, is in informal language in order to be put within 
the background limits already set previously. 

It was shown in (section 2. 21) that each plane (n ) determines the quadruple 
[a, 0, y, 6] up to a common constant of proportionality. That means that the 
planes [ a, 0, y, 6] and [Xa, X/3, Xy, X6], X/ 0, which are really parallel to each 
other in the Euclidean space, are considered one and the same plane in the 
projective space, or as one says in formal language, they belong to the same class*. 
Therefore, each one of the planes u 0 and ? r determines a class in the projective 
space, namely, they determine the classes FXa 0 , Xft>, Xy 0 , X6 0 ] and [P a , fJ-0, /J-y, P-6] 
respectively. Since the general Euclidean metric, where the metric d(7T 0 ,77 ) comes 
from, is a projective metric(see section 2,24), d(iro,ff) is a suitable metric for 
expressing the convergence of a sequence of "points" in the projective space, i. e., 
d(ir 0 ,ff ) can be used to express the convergence of a sequence of classes. 
Consequently, d(rr 0 ,n ) expresses the "distance" between the classes [X%, X0 o ,X%, 

X^] and F pa, pj3, py, p 16], independently** from the values of X and p. Thus, that 
"distance" is the same for any pair of planes from the two classes. Just for 
convenience one may choose the pair of planes for which 6 0 /ft = 6/p = 1. Note 
that the planes 6 0 = 0 and 6 = 0 are excluded from the respective classes. 

The above informal proof, indicates how analysis treats some intuitively 
appealing facts in Euclidean geometry. Another important thing here is the projective 
character of the Euclidean metric, 

*The "points" of the projective space are the classes [Xa, X0, Xy, X6], (X / 0). 

The class is completely determined if only one of its members is known. The 
others are obtained by varying X. 

**X and p keep their signs while they vary independently. 


33 





Consider now the two planes 

7T 0 : OioX! + j9bXg + y 0 Xa + 60 = 0, and (39) 

it : ax! + (SjXg + VxXg + 6 = 0, (40) 

as well as the lines 

l intersection of tt 0 and n 
£ perpendicular to l in tt 0 
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£ perpendicular to -t in it , and 
77 perpendicular to n. 


Without loss of generality, consider these lines being concurrent at the 


point P of l. 

Let the following notation be introduced 

x;. 

Xa, 

xS 

coordinates of P, 

V 

V 


direction-cosines of l, 

Y 

Y 

°« 

direction-cosines of £ . 

Y 

Y 

C c 

direction-cosines of £, and 

a , 

b , 

c 

direction-cosines of r). 

V 

V 

V 


Then the equations of and tt, as the planes which pass through P and are parallel 
to the pairs of lines [£, L] and [£, l~) respectively, are 


x x - x[ x a - x£ Xa - X| 



(41) 


and 


rr : 


xi -X{ 




xs-Xl 

Y 




= 0 . 


(42) 


t t l 

Comparing (39) with (41) and (40) with (42), the following is obtained: 


an 


& 


So 


70 = 30 


h c„ 


c 

4 4 

b„ c „ 

, a - s 

c c 

3. C. 

l l 



°c a s 

p 0 

, 0 = s 

C c "e 

p Q 

l l 



a_ b. 


& b 

4 4 

a. b a 

, y = s 

C 4 

a„ b. 

K/ ■L 


l l 


(43) 
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x; 

x p 3 

x£ 


x? 

x£ 

X§ 

$0 ~ _s 0 

a 4 

b ? 


, 6 -s 

a c 

b c 

C c 


a i 


c l 



\ 

c t 


wnere s 0 and s are constants of proportionality. Notice the presence of the point 
(Xj, X^, X$); Q is on t, generally different from P. The following theorem is 
useful at this point. 

Theorem . The square of the sine of the angle between two directed lines is equal 

to the sum of the squares of the two- rowed minors which can be formed from the 

matrix constituted by the two sets of direction cosines of the lines. 

If u and v are two lines with direction cosines [a , b , c ] and [a , b , c ] 

U H H V u v 

respectively, the above theorem proves that 


sin (fi, p) = sin*0 


a a _ 


b c 

3 

c a 

2 

a b 

4 (i 

+ 

4 4 

+ 

4 4 

b c 


c a 


a b 

V V 


V V 


V V 
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Proof . Consider the identity 

(a* + b® + cf)(af + b® + cf) - (a^a + b^a + c^g) 9 = 



b/ 

Cl 
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Cj a : 

3 

a i 

— . 

b 8 

eg 


Cg ag 

+ 

a 3 bg 


(45) 


Applying this to the direction cosines of the lines fiand v, the following is obtained. 


1 • 1 - (a a + b b 
4 v n v 


+ c c ) 8 = 
4 v 


b 

c 

3 

C 

a 

s 

a i bj 

4 

4 

+ 

4 

4 

+ 

b 

c 


c 

a 


a b 

V 

V 


V 

V 


V V 


But 


aa + bb + cc = cos (u, v) - cos 0 
H V fJL V fJL V ^ 

Then the left-hand side of the last relation is sin 3 0, and thus the theorem is proved. 
By virtue of this theorem, relations (43) yield 


0?a + ffi + y% -= s? sin 2 (4,t) = s 0 2 sin 3 90° - s£, 

and 

a* + 0® + y® = s 3 sin® (£, l) = s 9 sin 9 90° = s®. 
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Thus 


Pb = Sq and p = s. 


Relations (43) yield the normalized coordinates of the planes n 0 and it 
in terms of the direction cosines of three characteristic lines (t), (£), and (£). 
This representation is very useful for it offers the means to transfer from plane 
to line coordinates and vice-versa. 
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One may introduce the last relations into (38) to express d 9 (TT 0 ,7T ) in terms of the 
directional cosines of the lines (£), (£) and (£). However, the discussion continues 
at this point in plane coordinates. 

The expression (38) of d*(TT 0 ,7T) may be written 

a 0 a + W + y 0 y 




_ 2 


Pd 


= 1 + 1-2 


p 3 

a o a + p o 0+y o y 

PoP 


PoP 


or 




( 47 ) 
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As it is well-known, the vectors [a 0 ,£ 0 , y 0 ] and [oc, jS, y] are 
perpendicular to the planes n 0 and n , respectively. But each plane has 
two faces, and if [<*,$, 7] is pointing outward on one face of it, then the 



~ n [®0> Vo] > ~1To ~ 

HO Pc 

V 5 ~ n = ^ [-a, -/3, -y] 

Therefore, when one talks about "distance" of two planes, he actually 
means "distance" between two of the four faces. When the equation of a 
plane is written down, the face of the plane to be considered is defined. 
Consequently, that equation is for one of the faces of the plane. This 
remark sounds probably trivial, but it is of importance at this point to 
the advancement of the present analysis. 
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Let 0 be the angle* between the vectors [ ot 0 , &» 7o ] and ja,0,y], 
that is the angle between the faces corresponding to these vectors. Then 
one has 


“a. . + A . A. + la. . 2. 

ft) P ft P Po P 


cos 0 , 


and the relation (47) is written 

d 3 (TTo , TT) = 2(l-cos0) = 4 sin 2 (48) 

This is the "distance" between the faces [a 0 > A)> To] and [o:,/3,y]. 
Consider now the "distance" between the faces [ & , P, y] and [-Q! 0 , -ft,, -y 0 ]. 
Equation (47) gives 

,fi, _ (-ao)a > (-&>g > (-ro)r , 

cr (-ff 0 , if) - 2 *1- ftp ; • 

or 

d 2 (-ir 0 ,ir) = 2(1 + cos 0) = 4 cos 2 ■“ . (49) 

Thus, one observes that given a face tt = [a:,P, y] one may define with 
the two faces of another plane ir 0 two different "distances. " As a special 
case, one can consider the "distance" between the faces of a plane TTo; 
that is 


d (7T 0 • “TTo ) = 2 1 - 


Q! Q(~ a a ) + A ? (-ff o ) + y(r r-a.) = 2d + 1) = 4 

Poft) 


i.e 


9 


d 3 (ff 0 , ~ TTq) = 4 . 


(50) 


Recapitulating the above formulas about "distances" between the faces 
of two planes one has 

*One may find the notion "angle" out of place in this context. However, 
this discussion aims at revealing the geometrical meaning of the "distance" 
between two planes, and this mixture is unavoidable. 
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( 51 ) 


d(ff 0 , it) = d(-fr 0 , -it) = /2(1 - cos 6) = 2 | sin^ | 

d(ffo , IT) = d(7To , -IT ) = / 2(1 — cos 0 ) = 2 | cos~ I 

d(TT 0 , n 0 ) = d(7T, IT) = 0 
d(7To, -TTo) = d(JT , ~ IT) = 2 

In order to comprehend better the discussion following this point, 
a geometrical (or trigonometrical) interpretation of the formulas (51) is 
introduced here. Consider the cylinder with axis the line of intersection 

(£•) of the planes tt 0 and tt and cross-section the unit circle. The cylinder 

is considered of infinite length (Fig. 9). 






Take a cross-section of the cylinder (Fig, 10), This is a 
unit circle with the traces (AF) and(BC) of the planes tt and ff 0 , respec- 
tively. Using elementary trigonometry it is obtained that: 


0 

d(tr 0 , TT) = d(-7To , -it)' = 2 | sin— | = 4 • area (OCL) = (AC) 
d(-ff 0 , it) - d(ff 0 , -tt) = 2 I cos ^ I = 4 • area(OAM) = (AB) 


(52) 


Now, one has in his hands a very simple model in order to study the 
"distances" between the faces of two planes. The four faces +ir, -ff, +ff 0 , 
-ffo are modeled by the points A, F, C, B, respectively on the circum- 
ference of the unit circle, and the "distance" between any two of the faces 
is modeled by the chord joining the corresponding points. One may check 
all the formulas (51) through this model. Thus 


dfro, tt) = (AC) 


d(-ffo.-ff) = (BF) 
d(-ff 0 , ff) = (AB) 

d(-ff 0 , tt) = d(7T 0 , -IT) = 2 | cos | | , 

d(7T 0 , -IT) = (CF) 


d(7T 0 , TT) = d(-ffo, -ff) = 2 | sin | | 


d(ff 0 , -ffo) = (BC) = 2 
d(ff, -it) = (AF) = 2 


^>d(1T 0 , -ff 0 ) 


d(ff, - tt) - 2, 


and obviously 


d(ff 0 , Tb ) = d(ff , n) = 0 
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Remark . One thing that should be stressed here is the linear representation of the 
"distances" between the faces of planes. It should have been noticed that it is not 
necessary to inverse the circular functions involved, and express the "distances" 
in terms of arc. 

The two different "distances" which can be defined between the faces of two 
planes are not independent from each other. They are related through the relation 

d 3 (tt 0 , TT) > + d 2 (-7r 0 , 7) = (AC) 2 +(AB) S = 4 ( 53 ) 

Having the distances between faces of two planes, naturally springs up the question 
as to whether there exists a way to treat each plane as a unit and find a unique 
"distance" between two such units. That "distance" has to be a Euclidean metric 
and consequently a projective metric (see section 2 . 24 ). The above-employed 
representation for the "distances" between faces of two planes (Fig. 10 ) shows 
so clearly the existence and uniqueness in magnitude of such a "distance", that it 
hardly needs any further explanation. Thus the "distance" between the planes 77 and 
n 0 is given by the distance (AN) of A (face +77) from the hypotenuse (BC), (trace of 
the plane on the plane of the cross-section). As one checks immediately the 
distance of A from (BC) is equal to the distance* of F (face ~n) from the same line 
(BC), as well as with the distances of B (face -n o ) and C (face + tj 0 ) from the line (AF), 
(trace of the plane 77 on the plane of the cross-section). Therefore, the "distance" 
between i t 0 and 77 is uniquely defined. Let it be called D(n 0 , n)- It remains to be 
proved that D(t7 o , 77) is an Euclidean metric. To this end, the geometrical meaning 
and the connection of D(tt o , 77) with the "distances" , 77) and d(_T 7 0 , 77) will be found 
as intermediate results. From Fig. 10 one obtains 

area (ABC) = & (AC) (AB) = | (BC) (AN) . 

= | . 2 . D(t7 0 , 77) = D(t7 o , 77) . 
or, by recalling ( 52 ) and ( 51 ) 


*What has been said for A could be repeated for D. Because of the antisymmetry of 
the cross-section the discussion is given on the one half of the figure only. 
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D(tt 0 , tt) = \ (2 sin | ) (2 cos | ) = sin e . 

Thus 

D( 7 T 0 , n) = sine , (54) 

and also 

D(n 0 , n) = area (ABC) = § d( Uo , n ) ■ d(_rr 0 , it) • (55) 

The extreme values of D(rr 0 , tt) will be given here without any recourse to the 
properties of sine. From (53) one concludes immediately that since the sum of 
d 3 (n 0 ,r r) and d 2 (_ 7 r 0 ,ri) equals constantly to 4, their product becomes maximum 
when d 2 (n 0 , tt) ~ d 2 (~tt o ,tt) = 2. Thus the maximum of D 2 (tt c » tt) = ^ d 2 (t r o , tt) • d 2 
(-ir 0 ,tt) is 

max [D 2 (tt o ,tt)] = 1 

Since D(t r 0 , n) is taken positive by definition, it can be written 

max [D(tt o ,7t)] = 1 or D 0r o .7r)^l> ( 56 ) 

that verifies the known result jsine <1. One can verify all the known properties 
of sine- 

In order to be shown that D( 7 r 0 7 r) is an Euclidean metric it is sufficient to recall 
the relation (55) , i.e., 

D 3 (rr o . n) = ^ d 2 ( TT 0 , tt) ■ d a (-TT a , n) 

where cP(Tr 0 ,t t) and d a (_Tr 0 ,tt) are Euclidean metrics, as it is known from the 
preceding discussion. But it is worthwhile to put D 2 (tt 0 ,n) in the general form of an 
Euclidean metric ( see section 2.24, relation (18) ). One may proceed as follows: 

D 3 (tt 0 ,t t) =4 d 2 (tt 0 , tt) . d 2 (~tt 0 , tt) =~ 2 (l-cose) 2 ( 1 +cose) 


= l-cos 2 e 



PoP 
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The last expression may be written 


D 3 (tt 0 ,tt) = l^U 


( 57 ) 


(58) 


where 
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(59) 


This shows that D 3 (tt 0 , ir) is a special case of the general egression (18) for a 
Euclidean metric between points (see relation (32)). Thus, the metric D( 7 t 0 ,tt) 
will be called dual metric of the usual Euclidean metric. As it was proved, 

D(t7 0 , rr) renders the shortest "distance" between the planes n 0 and n (considered 
as units), in the same way as d(P 0 , P) renders the shortest distance between 
the points P 0 and P. 

The above expression (57) of D 3 (tt o , 77 ), is in terms of plane coordinates 
actually normalized plane coordinates. Two more analytic expressions of 
D 8 (n 0 tt) very useful for the adjustment problem will be given here. One is in 
line coordinates (directional cosines) while the other is in line and plane coordinates 
For the first, consider the lines (|) of rr 0 and (|) of the n perpendicular to the 
intersection line ( 1 ) of n 0 and 77 . Then relation (44) yields: 


D 3 (rr 0 , tt) * sirfe = 

b s % 

2 

+ 


2 

+ 





c ^> 




For the second, consider again the line ( 5 ) of n 0 (the same as above), and the 
vector [a, /3, y], (which is perpendicular to tt). Then 
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Recapitulating the results yielded by the above discussion about the 


"distance" between two planes one has the following: 


« I Qri/s i = (n, tt)CE \l 


(i) Given the equation otx^ + /3^ + -yXg + 6X4 = 0 (of< aplaheiir ,(in)homogeheojxs v G \S 
coordinates) or class [tt] of planes is, determined, whose! representation vector 
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is |Xa, X/3, X-y, X6], where X^ 0-1 F ; orX>0 the^subclass [-f-TT^is jobtained 

I B ss | j ® Y 1 I Y 

determined by one of the faces of_the„plane TT;.namely, ^theJace_+n_=_[d,’ 
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For X<0 the subclass [-irj is obtained whichUs' determined'by^theJoiher face of it, 
namely, the face -77 = [~oc, -/?> ~yj-: 
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d(_rr 0 , +tt) = d (+tt 0 , _tt) = /2(l+cose) - 2 j cos'g 
d(+TT 0 , -tt 0 ) = d(+TT, -tt) = 2 
d(+TT 0 ,+TT 0 ) = d(+TT,+TT) = 0 
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These are Euclidean metrics, and express "distances" between faces of two planes, 
(iii) There exists a uniquely defined Euclidean metric, D(tt 0 , tt), which expresses 
the shortest "distance" between two planes tt 0 and tt. The following expressions for 
this metric are equivalent: 


1/ D(rr 0 , tt) 
2/ D(ri 0 i tt) 


3/ D(tt 0 , tt) - 

4/ D(rr 0 , tt) 
5/ D(tt 0 ,tt) 


sine I , 
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where [a^, b^, c^] and [a^, b^, c^] are the direction cosines of the lines (f) and (£), 
which lie in the planes rr 0 and tt respectively and are perpendicular to the intersection 
line (t). 

Among all these expressions the last one is of particular interest for the 
adjustment problem (see section 2. 22), and it will turn out to be very convenient. As a 
matter of fact, this is a linear expression with respect to a^, b^, c^. and J2f, $ , , 

and fits the adjustment problem as it was formulated in section 2. 22. Thus the 
required expression for the regression model has been found. The following subsection 
is devoted to the statistical analysis of this model. 


2.26 Statistical Analysis of the Adjustment Problem , 

mgtTTTTZ aiij 1 ~ ; S "LV — — J lUl SIl T .SQXJlI flOf J083185Xn . 

The geometrical analysis of the optical observations mode (see section 2. 1) 

sril SB 80X8 a 3 ill 8l JbiIj lO ,381111 XJOllOea'ISlill Qxil 10 JOB Sill OJ Qfixl 

revealed the type of the statistical model, namely, a linear regression model. 
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Although the enunciation of the adjustment problem (i. e. , what is the least- 
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problem, the analytical expression of the model was not so obvious. Being 

.bsaailsa ai nolislrxo muaiilqo asiBijpa IbbsI gniwollol 9xii dob 
convinced of the linearity of the problem, one was expecting^ linear analytical 

expression. Duality principle offe^^^qipjcan^ifeo fjrytj ?®t the suitable analytical 

expression for the linear regression model. The fact that dife knows the statistical 

model does^not discharge him from being obligated to discuss it stochastically 

and to appraise the vaTidt ty-o f The s ta ti s tic a 1 assumptions which accompany the 

| — - — (V \ 

model, under the conditions of the problem. Thrs-suhsection-aims at a general 
statisticaLt]|eatment/Of the problem. 

VI V \ 

2. 261 Statistical Model - — 

Two| ground stations Pi and P a ar^~ ca-rxyipg^^ s^mudtaneou s optical 
observations of a set~of “ sa'f elRfe - positron - poinFs~Q~V Q p . The 

observed Ouantdtie&sb^^ Pi at the event [E 3 ,Qj | tj ](qee section 1) are 


equivalent \|o the direction cosines^f^p^^^^j^i j ) of the direction PiQj. 

At each event]Ltwo m cbx^ctiojis.^EiQ.j^Jid-Tl a l3-j m ^^5^^Bnaned . which 

"0“- p<[ 6+ C x jV + r x.a-f,x .Q : Jt \ 

are not independent from each otner, for they meet at Qj. Th^se two 
directions intersecting at Qj, determfneSh^ plane rr i which is defined as 
passing from Q 3 and being parallel to the directions PiQj and P a Qj(see section 
2.1). Thus, there 1 is-*h?)s£t ,(p) T inde^endgnt^ v/ 

determined from each other. One might think that since the planes tt^ are 
passing through the line (P a P 3 ^,i&hpy^cofold bje ..cSnsic^e^ed ap^-piembers of 
a pencil of planes with axis, The £ °line (PiP 2 jt ^t such a thing would 
kill the problem statistically, for the places so -c|ef inex|, , are^, not inde- 
pendent from each other any more. Further discussion on this is post- 
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poned until later. These planes determine Cs p = p(p - l)/2 independent 
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intersection lines. The line (P v P a ) is considered as the fitting straight 
line to the set of the intersection lines, or that is the same as the 
least -squares straight line through the set of planes ttj (j = 1, 2, 

Let (4) be the line (P x P 3 ) with direction cosines [a^,b^,c^]. Then 
according to the discussion in section 2. 22, it corresponds to the line for 
which the following least squares optimum criterion is satisfied. 


V 

^ dP(tt oJ , rij) = minimum. 


(62) 


J=i 



where D(n 0j ,TT 5 ) (see section 2. 25) is given by (61), i. e. , 

y£v 

Pi = [Off + 0j 3 + yf . 

One may write 

Dj = D(£, rr J ) = D(tt 0 j, tt j) = 0f N j a^ + &,}bg + y N} c^ , (63) 
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where 


a , „ §_i_ y, c _ 6, 

s "‘ • 


are the normalized coordinates of n i , 
If (63) is written in the form 


(64) 


D 


= aN * a g * ^ Njb £ + yN J C g + 
[a> +b t 3 + c *] 2 


where as is known 


[a> + b t 2 + c t 2 ] = 1 


( j = 1, . . . , p) 


(65) 


(66) 


then one recognized immediately the formula for the distance of a point 
with coordinates [ ^ N j, #n j,y N j] from the plane, the perpendicular vector 
of which is [a^, b^.c^]. This is very remarkable because the problem 
in question is thus reduced to a very well-known case, namely: find the 
least-squares plane through a set of points. Hamilton's presentation 
of the least squares solution to this problem will be adopted here [18, 19]. 


Let 

and 


m = [a., b^, c f ], 

X J - [®N J > 0N j , yN j] * (j — 1» 2, . . . , p) 


(67) 

( 68 ) 


be the coordinate vectors of the "plane" to be fitted and a typical 
observed "point"* respectively. Then (65) and (66) may be written 


Dj = 


mX) + o 
[m m T ]2 


m m T : 


mxj, (j = 1, . . . , p) , 


(69) 

(70) 


*The quotation marks are justified by the fact that the notions point and 
plane are symbolic there. 


49 



The equations (69) are the so-called observation equations while (70) is 
an absolute constraint. As far as statistical assumptions go, it is 
assumed that j8 fJ j and yNj are independent normal random variables. 

Before proceeding to the least squares solution, some comments on the 
statistical model are set forth. 

2.262 Comments on the Statistical Modal 

It certainly has been noticed that the "observed quantities" in the 
above adjustment are the coordinates a NJ , fin], y N j of the plane tt i . 

One might wonder how far away these quantities are from the quantities 
which are actually recorded by the observer. In Principle, the quantity 
which enters the adjustment as observed, represents the observational 
unit , so to say, that is the event which is repeated during the obser- 
vations and which may be considered independent from the preceding and 
following one. The question now is, what is the observational unit in 
the optical observations mode? The "instrument" of these observations 
is not each one of the two instruments at the observing stations, but 
both of them, and the "observed quantity" is not each one of the two 
directions, but both of them together as a unit. That unit is nothing 
else but the plane defined by these two directions which, by fact, meet 
at the satellite position point. Obviously, this is the observational 
unit, or the event, if one wishes, which is repeated during the observa- 
tions and may be considered independent from the other such events. 

Therefore, the quantities which enter the adjustment as observed 
must represent the planes which are defined independently from one 
another and each by the corresponding simultaneous directions. Thus, 
the above adjustment cannot be simplified further without violating basic 
statistical principles about adjustments. But there are many ways to 
represent one and the same plane. This matter was discussed in section 
2.1, where it was concluded that the only statistically permissible 
representation of the plane tt 3 (see Fig. 12) is from the two 
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Fig. 12 


intersecting directions PiQ 3 and P 2 Q 3 . If one considers the planes tt 3 , 
forming a pencil of planes with axis the line (Pi P a ) v then the independence 
of these planes disappears, and a mathematical model for which this is 
the case, does not yield an adjustment of this problem. From the above 
discussion it became clear that there was no other way to look at this 
problem, (i.e. , two unknown ground stations observing simultaneously 
directions to a set of satellite position points). That means one must 
look at (PiP 2 ) as the fitting line of a set of lines (all possible inter- 
sections of the planes tt 3 ). i n the following of the present commentary, 
a mathematical model for the same problem different from the above which 
appeared in geodetic publications will be discussed briefly [2]. The 
mathematical model given here there (somehow arbitrarily) is the following 
(see Fig. 13). 

Direction PiQ 3 : 

x 3 ~ Xi - r 13 a X3 = 0 
y 3 “ Yi " r 13 b 13 = 0 
z l " zi - r uCi 3 = 0 
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Direction, P s Q i: 

x j - Xz - r 3j a2j =. 0 
Y J - Ya - r 2J b 3 j = 0 
Z J ~ Za. r sj c 3J = 0 



Fig. 13 


There is no question that they are necessary, but as it turns 

they are not sufficient for this problem. One may observe the following 
points. ^ 

After elimination of the satellite coordinates, the above conditions 
are reduced to tllree ^ 

f h 6 three ^“Itent equations have the form 

- xs + rijaij - r sj a 2 j = 0 

(only the x-component is written for brevity). 

in th° ne h migl r tWnk ‘ he PreSeaCe °' ‘ he COOrdiK ‘ tes *nd ‘he distances 

o, 11 Td ‘ S '" tri " Sl °- ^ ‘ hiS iS n0t thC — • matter 

feet, the distances r u and r a] are unknown and they become known 




after the introduction of the length ri a , through the relations 


and 


Then 


or 

(sin0 3 ) aia +(sin9 a ) aij - (sin0i) a a3 = 0. 

The last one is independent from the coordinates of the points Pi and P a . 

' 

Another thing which must be pointed out is that the directions PlQj 

and P a Qj are introduced independently, while they are not independent for 

they intersect at Qj . This dependence is expressed by the equation of 

— ♦ — * 

the plane through Qj and parallel to PiQj and P a Qj, which is not 

mentioned anywhere there. These two points show that the above model 

does not give the required adjustment. It is not easy to visualize what it 

— * “ ♦ — ■* 

does. It seems that after such an adjustment each triple [PiP a , PiQj» PsQj] 
is parallel to a plane, but there is no connection between these planes. 

That model cannot be identified with any of the known statistical models. 

It belongs to the category of the so-called generalized least squares. 

There appears to be no objection against the generalized least squares, 
provided the employed conditions are proved to be not only necessary but 
sufficient too; sufficient in the sense that they define the configuration 
uniquely. Unless the sufficiency of the conditions is proved, the results 
are very precarious. The fact that there is no rigorous statistical 
foundation for the generalized least squares is an additional reason for 
some people to keep their reservations with this method. In the last 
analysis, why should one have recourse to an unconvincing generalized 
least squares solution, while he has a clearer, fully convincing solution 
of the adjustment problem ? 

53 N 


r U = 


sin 0 c 


r a j = ria 


sin 03 

sin 9 t 
sin 03 


XL-USa + 8 in, 9 a a _ 8ip ,§l 
r ia sin 03 13 sin0 3 
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2.263 Least-Squares Solution ; 

Equations (69) and (70) constitute the mathematical model of the problem. 

The least-squares solution to this problem is very known. As a matter of fact 
there are more than one approaches to the solution of this problem [17, 18]. 

Since this mathematical model is very common in statistics there exists a 
vast bibliography about its statistical analysis (computing techniques for estimation, 
hypothesis testing, etc.) which is very beneficial for the geodetic problem. Once 
the problem of the optical observations mode is reduced to that one above, the task 
of the present work in this direction has been accomplished. Nevertheless for com- 
pletion purposes, Hamilton's method of solution will be given below. 

Hamilton's solution to the problem of fitting a "plane" through a set of "points" . 

The following matrix notation is introduced 


d = (D lt n,, • • 

• , D p ) , 

(lxp) 


(71) 

X = (x x ,x a , . . 

• • ) . 

(3xp) 


(72) 


Mjj. ( 3p x 3p) variance-covariance matrix of X 
(pxp) variance-covariance matrix of D 

W = M d -1 weight matrix of D [to be discussed later] 

Point estimation. 

Proceeding formally one obtains the function to be minimized, i.e., 

§= D W D T - k(mm T - 1), (73) 

where k is a Lagrange multiplier for the constraint mm T = 1. From (69) and 
(70) one has 

D = (mx^, mx a , . . . ,mx p ) = m X (74) 
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Then (73) may be written 

$ = m(X WX T ) m T - k(mm T - 1). 


The least-squares estimate of the unknown vector m is given as it is known by 
the solution of the following system : 

= 0 

am 

or 

(X WX T -kI 3 )m T = (C-kI 3 )m = 0 (75) 

But the last equation shows that m T is an eigenvector of C corresponding to the 
eigenvalue k. Since C is a symmetric positive definite matrix it has three real 
eigenvalues, and let them be in the following order 

^ i ^ ^ ^3 • 

Then in T corresponds to k = min \-X l . By some standard numerical 
procedure one may obtain the minimum eigenvalue. The eigenvalues of C -1 are 
the inverses of those of C, and the maximum eigenvalue of C _1 corresponds to 
the minimum of C. Since the numerical procedure to find the eigenvalues starts 
with the maximum eigenvalue one is looking for the maximum eigenvalue of C -1 , 
whose the inverse is the minimum of C. Given the fact that C = X WX T is a 
(3x3) matrix there is not much numerical trouble involved in these manipulations 


Nature of the Weight Matrix W. 

In this problem the variance covariance matrix of X, i. e. may be 
considered diagonal with (3x3) diagonal blocks. That because, the planes iTj , 
(j =■ 1,2, . . . , p) may be considered independent. Nevertheless in this general 
treatment, a full variance-covariance matrix is taken and let its typical 
element be denoted by IV^ 13 , i.e., 



. . M x ’ P 

— P » 1 

— p ,p 

M 

. . M ’ 

Y 


M*’ 3 is (3 x 3) 


(76) 
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Equation (74) may be written 



Then the formula for the propagation of error yields 

w" 1 = M d ’a -30 M x “M T 


( 77 ) 


(78) 


The typical element of M d is 

— 1 J _ lj _ T 

M d = m M x m (79) 



or 



Since the weight matrix W is a function of m the equation (75) is non-linear in m. 
Therefore an iterative procedure is followed. One assumes a vector m, calculates 
the weights from (78), derives a new value for m from (75) and then recycles. In 
cases of two or three dimensions, like the present one, the process converges in 
no more than two iterations. 
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Remark . Because the weights depend on m and an iterative procedure is 
used it is important to note that the minimum eigenvalue of C does not necessarily 
correspond to the best solution. When the solutions to the eigenvalue equation (75) 
have been found, the quantity D W D T should be computed for each eigenvalue. 

If the ellipsoids of error for the individual points are very eccentric, it frequently 
happens that the minimum D WD ? does not correspond to the minimum eigenvalue [20]. 
This is particularly likely to happen if the initial choice of m for the iteration is 
far from the correct one. 

Covariance matrix of m ( Mm ) . 

From (74) one obtains 

D T = X T fn T or X W D T = X WX T m T = C m T 
and 

m T = (C -1 X W) D T (82) 

By the formula for the propagation of error the last relation yields 

Mm = (C _1 X W) M d (WX T C -1 ) 

= C _1 X W M d M d _1 x\x WX T ) 1 

= C _1 (X WX T ) (X W X V 1 

or 

Mm = C“* (83) 

Testing hypothesis ■ 

Consider the hypothesis H: m o = m. Assuming that the distribution of D 
is normal then the quadratic form 

Q 2 = (m 0 - m) C(m 0 - m) T (84) 
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is distributed as 3F 3tP , where F is the well-known variance-ratio distribution 
and 3 is the number of unknown parameters. If Q? is less than the tabular value, 
the hypothesis is accepted. That means may characterize the "plane, " which 
is the line (£) in this problem. Since 


lim 3F 3 , p = Xa 

P~>CO 


for large p one may use \ a test instead of F. 

In every statistical problem, after the statistical model has been chosen, 
the goodness of the estimation and the statistical inferences depends upon the 
weights assigned to the observed quantities. For the present problem, the 
difficulty is about assigning the weight to the observation vector"xj. By recalling 
equations (2) - (5) (see section 2.1), one has 


or 


X J = 

[o'nj . 


- i fc 

» Pi ’ Tj ] 

X J = 


b u c u 

> 

c u a \j 

y 



^3] C 2i 


C 2J a Sj 







1 



Pi = 

[o' 

f + P? 

+ yf 

f • 




(85) 


where 


It is interesting to note that p 3 is the sine oF the angle which is formed by 
the two observed directions at the satellite position point (see relation (44)). 
If the observed angular quantities are the declination (6 ) and the hour angle (h) 
the following relations hold 


•i j 

=. cos6 1J 

cos h^ j 

^ i 

= cos6 2;) 

COShg 

1 j : 

= cos ] 

sin hj j 

^ j 

= COS 6 2 J 

sinhg 

I j 

= sin 6 XJ 


°Si 

= sin 6s j 



( 86 ) 
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Then 


b u c i 3 cos 6 X j * sinh 1 j , sin6 XJ 

ot } = = 

b 2 j C^J cos6 SJ * sinhgj , sin 6^ 

a u C XJ cos6 u * cosh 1Jf sinfi^ 

~~ = ~ * 

32J c 2J cos 6 S j * coshg j , sin 6 a 3 

a i j b XJ cos 6 x j cosh l3 cose^ sinh xj 

Tj = 

b 2 j cos 6a } cos hg j cos 6 a i sinh^ 


(87) 


(88) 


(89) 


Pj = sin(P 1 Q J , PgQj ) = sin6j (cos 1 0^ ) 

= t 1 “( a u ^ J +b iAj + c lJ c SJ ) a ]^ 
or 

Pj 2 = l-[sin6 13 sin6 SJ +cos6 XJ cos6 aj cos(h x j -hg } )] s (90) 

The above formulas are not too complicated to compute the errors of 

a N j , /3 N j , y NJ , by using the formula for propagation of error on the linearized 

expressions of the above when the errors of 6 1 } , h xj , 6 a 3 , hg 3 are given. It 

is stressed again tnat the observational unit of this problem is the plane tTj , 

— » — > 

and not the directions P x Qj and P 2 Q, individually. Hence, the coordinates of 
this plane are the "observed quantities" for the adjustment problem and these 
are those which must be assigned a weight. As it was inferred above, the plane 
TTj is the simplest observational unit and any mathematical model whose observa- 
tional unit is simpler than that plane does not adjust the observations of this 
problem. 
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2. 3 Determination of the Cartesian Coordinates of the Stations 
After the above adjustment one has at its disposition the direction cosines 
of the line (£) = (PiPg). In order to compute the coordinates, one needs the 
coordinates of either one of the two points and the length. Let x^yi.Zi be the 
coordinates of P L , and S 12 be the length of the distance (P^). Then 

Xg — a^ Su2 + X; 

y 3 - b^S l3 + >T (91) 

z 3 " c £ ^13 + z 1 

3. THE CASE OF MORE THAN TWO STATIONS 
Let P lf . P 8 , .... P t be a set of ground stations. Consider all the pairs of 
stations from this set, each one of which has observed simultaneously at least one 
satellite position. It does not matter whether some of the satellite position 
points take part in more than one pair. Let these pairs be classified 
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as follows: 


A 3 set of station pairs each of which has observed more than two satellite 
position points, 

A 2 set of station pairs each of which has observed just two satellite 

position points, 

A x set of station pairs each of which has observed just one satellite 
position point. 

Remark . It is possible, of course, one and the same station, say P k , to be 
found in all these three sets of pairs, but never as it is understood, two and 
the same stations can be found in two different sets of pairs. 

1st Estimation. 

For each station-pair of the set A 3 , an adjustment is carried out 
which yields the directional cosines of the line of that pair. For each station 
pair of the set A 2 , there is no adjustment problem inasmuch as there exists 
only one intersection line, which is taken as the line of that station-pair. If, 
for example, (P^Pa) is a station-pair which has observed two satellite points 
Qi and Q 2 (see Figure 14), the only intersection line is that of the planes 
TT\ and tt 3 (see Section 2.1). For the station- pairs of the set Ax there is no 
information during the first estimation. Thus, the first estimation consists 
of the direct information provided to the ground stations from the satellite 
observations, that is, the directional cosines of the lines of the station-pairs 
in A 3 and A 2 . 

After the first estimation one has additional information for the deter- 
mination of the directional cosines of the lines of the station- pairs, which 
may be called indirect observations. That is, in a second estimation of the 
direction cosines of the line of a station-pair (Px,P 2 ) one may consider not 
only the planes tt } corresponding to the satellite position points Qj observed 
by the stations Pi and P 2 , but also the planes II^ = (PxP 2 Pi)(i = 3, 4, ...), 
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provided the directions Pi Pi and Pa Pi have been estimated in the first 
estimation. Consider, for example, the simple configuration of three 
stations P^Ps.Pa (see Figure 14), where the lines (PiP 2 ), (P 2 P 3 ) and 

(P 3 P 1 ) have known directional cosines from the first estimation. The 

-» -> ' -* 

directions P 1 P 3 , P 2 P®, PgPi as they result from the first estimation, 
are not in general coplanar. One might consider here the coplanarity 
condition of three directions, i.e. , 


ai a 

bi 2 

Cl 2 



a 2 3 

ba3 

C23 

= 0 , 

(92) 

agi 

boi 

C 31 




and adjust the results of the first estimation according to this condition. 

This way is followed in [1].. But taking the first estimation results and 
enforcing them to verify condition (92) it is like killing the problem statis- 
tically, and abandoning the principles adopted in this treatment for the case 
of two stations by implicating statistical models for which theories about 
statistical inferences have dubious validity; the function F (aij , b t j , c tl ), 

(i, j = 1, 2, 3) has nine variables and it is of third degree. Therefore, 
by principal, the adjustment here will not proceed this way. The natural 
way to proceed is that of a second estimation, as already quoted above. 

2nd Estimation. 

The second estimation is similar to the first, and only the number of 
observations has been increased, by taking into account not only the direct 
observations (which yielded the first estimation results), but the indirect 
observations too. Consider a station-pair, say (P x , P 2 ), and let the following 
sets be defined 

B — [ 7T ] — [ TTi , ...,7Tj,...] 

c = [ii] = [ii 12 . . . , n 13 , . . . j , 

where B is the set of planes corresponding to the satellite points Qj having 



been observed by the station-pair (Pi, P 2 ) and C is the set of planes 
_ 

Hi s (PiP 2 Pi) corresponding to the ground stations Pi, such that Pi Pi 
and P 2 Pi have been determined in the first estimation. It is observed 
that in the second estimation one obtains the direction cosines of a 
station- pair whenever the set C is not empty, (the set B is not empty 
by fact). Thus, while one could not estimate the direction cosines of a 
station-pair in the set Ai in the first estimation, he might be able to do 
that in the second one, provided the set C of the station-pair in Ai is not 
empty. 

Note. One should be aware of the definition of the plane n i2 . It is defined 
as the plane which passes through the station Pi (upper index) and it is 
parallel to the directions Pi Pi and P 2 P . 

After the second estimation a third one should follow in which the 
observations consists of the direct (original) and the indirect ones resulted 
from the second estimation. Thus, an iterative procedure is established 
which obviously leads to the coplanarity condition. The iteration steps 
continue until no significant difference exists between the results of the last 
step and the one immediately preceding it. 

After that the adjustment is complete and for the Cartesian coordinates 
of the stations one has to introduce the coordinates of one station and the 
length of the distance between two of the stations whose direction has been 
determined during the adjustment. 
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